Effect of spin on electron motion in a random magnetic field 
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We consider properties of a two-dimensional electron system in a random magnetic field. It is 
assumed that the magnetic field not only influences orbital electron motion but also acts on the 
electron spin. For calculations, we suggest the trick of replacing the initial Hamiltonian by a Dirac 
Hamiltonian. This allows us to do easily a perturbation theory and derive a supermatrix a model, 
which takes a form of the conventional a model with the unitary symmetry. Using this a model we 
calculate several correlation functions including a spin-spin correlation function. As compared to 
the model without spin, we get different expressions for the single-particle lifetime and the transport 
time. The diffusion constant turns out to be 2 times smaller than the one for spinless particles. 
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I. INTRODUCTION 

The problem of electron motion in a random magnetic 
field (RMF) has attracted considerable interest in recent 
years. One of the main questions is how the localization 
scenario is modified in comparison with the usual disor- 
der problem. It is also thought to be relevant for the 
behavior of composite fermions near half-filling in the 
fractional quantum Hall effect .El In a recent experiment 
for a RMF system a magnetoresistance similar to that of 
the quantum Hall system was observedfl 

Analytically, this problem- has been discussed us- 
ing bothj-thc diagrammatical and the supersymmetry 
method.aOu It turns out that the system belongs to the 
usual unitary class, which leads to localization in two di- 
mensions (2D) unless the random magnetic fields B cor- 
relate over a very long distance (1/q 2 dependence for the 
(B q B- q ) correlation). As concerns a short range correla- 
tions of the RMF, one can derive a conventional nonlinear 
a model using a standard procedure (see, e.g., Ref. ||). 
With long range correlations one can derive first a bal- 
listic cr-modelfl Integrating out nonzerQ jiiarmonics one 
comes again to the diffusive a modeLU'EKI Only if the 
correlation of the magnetic fields obeys the 1/q 2 depen- 
dence, one may get something different (antilocalizationA 
because, in this case, a new term in the a model appears.u 

A standard Hamiltonian used for the RMF problem 
has the form 



H 



1 / e \ 2 
n —[P- -A(r) 

2m V c w 



(1) 



where e and m are the electron charge and mass, and 
the vector potential A(r) corresponding to the magnetic 
field B (r) should be averaged with some weight. 

Eq.(yJ) describes electron motion in a magnetic field 
neglecting interaction of the magnetic field with the elec- 
tron spin. This is a good approximation for GaAs het- 
erostructures where the Zeeman splitting is very small. 
In models of composite fcrmions spin is absent at all 



and therefore the Hamiltonian Ho, Eq.([y), is sufficient 
for proper description. 

Nevertheless, the question about the character of elec- 
tron motion in a RMF acting also on the electron spin 
may be interesting on its own. Generally, the interac- 
tion of the magnetic field with spin is not smaller than 
interaction with the orbital motion. Moreover, in a ho- 
mogeneous magnetic field the Zeeman splitting for free 
electrons is equal to the distance between Landau levels 
(see, e.g., Ref. ||) and one may ask if this degeneracy may 
show up in the inhomogeneous field. To the best of our 
knowledge the problem of electron motion in a random 
magnetic field with both orbital and spin interaction with 
the magnetic field has not been addressed yet. 

In this paper we address this problem starting with the 
Hamiltonian H 



H = Ho- ^BtrB(r) 



(2) 



for g — 2. Here, [Lb = e/2mc is the Bohr magne- 
ton, a stands for the components of Pauli matrices and 
B(r) = V x A(r). It is clear that one may not consider 
the orbital and spin interactions separately. This is in 
contrast to models describing electron motion in a mag- 
netic field in the presence of magnetic impurities. The 
effect of the Zeeman term has been examined in Refs. 

for the case of a scalar random potential and a small 
constant magnetic field and corrections to the conductiv- 
ity have been calculated. 

One should notice also that, while the interaction with 
the orbital motion is described by the vector potential 
A(r), the interaction with the spin is determined by the 
magnetic field B(r) itself. Therefore, averaging over the 
magnetic field cannot be a trivial procedure and devel- 
oping a proper calculational scheme may be interesting 
from the technical point of view. 

We suggest a scheme that has not been used in the con- 
text of disordered metals. Our idea to discuss the RMF 
problem for the Hamiltonian with the Zeeman term is to 
consider a more general Dirac Hamiltonian as a starting 



1 



point of the analysis. There are two advantages to use 
the Dirac Hamiltonian. First, the square of the Dirac 
Hamiltonian gives Eq.(||). Thus, we can naturally take 
spin effect into consideration for the analysis. Second, 
the Dirac Hamiltonian contains only the vector potential 
A(r) but not the magnetic field B(r). The dependence 
on the gauge field A(r) is linear, which simplifies the 
averaging procedure. 

The Dirac Hamiltonian has been used for the prob- 
lem of random Dirac fermions. This problem may be 
relevant rior degenerate semiconductors JIIL quantum Hall 
systemsjla and d-wave superconductors ,li3 and has been 
under intensive study. For these problems, the chiral 
symmetry of the Hamiltonian, which means the energy 
eigenvalues are symmetric around the zero point, plays 
an important role. In contrast, the chiral symmetry is 
not important in our analysis since we consider energies 
in the vicinity of the Fermi level and therefore far from 
the zero energy. 

We show below that this model can be mapped onto 
the conventional nonlinear a model with the unitary sym- 
metry. As the chiral symmetry is not important, this is a 
natural result from the viewpoint of the symmetry con- 
siderations. At the same time, the interaction with the 
spin changes the classical diffusion coefficient D entering 
the a model. We obtain somewhat different expressions 
for the single-particle lifetime and the transport time as 
compared to the RMF model without spin. The spin 
degrees of freedom do not change the conventional form 
of the current and density correlation functions. As a 
spin dependent quantity, a spin correlation function is 
calculated. 



II. SELF-CONSISTENT BORN 
APPROXIMATION 



a (5-correlated magnetic field. Nonzero values of k de- 
scribe a screening and we keep below an arbitrary value 
of K. 

Our main idea for considering the spin effects in the 
random magnetic field problem is to represent Green 
functions G 



R.A 



G 



R.A 



1 



E- H±iS 



(7) 



of the Hamiltonian H, Eqs. (hi) and (B), in terms of the 



Green functions 



R,A 

9 k 



R,A 

9k 



of the Dirac Hamiltonian n 



1 



k-U±i5 



Using the relation 



n R,A _ 

up — 



m ( 

¥ V 



R,A 

9 k 



A,R 
1-k 



(8) 



(9) 



where k — (2mE) 1 ^ 2 is actually the Fermi momentum, 
we achieve this goal. 

The Dirac Hamiltonian n is linear in the gauge field 
A(r) and the ensemble averaging can be performed easily. 
On the other hand, the usual Hamiltonian (|l]) includes 
the square of the gauge field, which makes the averaging 
procedure more difficult. At the same time, just neglect- 
ing the term A 2 may be dangerous because this violates 
the gauge invariance. In our case, we can keep the gauge 
invariance at any step of calculations. 

In order to demonstrate how our scheme works we 
calculate first the average Green functions {g^.' ,A ^ us- 
ing the well known self-consistent Born approximatiop. 
Summing only ladder diagrams in the standard waytil 
we obtain for the Green function 



Before deriving the a model let us demonstrate how 
one can calculate the one-particle and transport lifetimes 
using the Dirac representation for the Hamiltonian H, 
Eq.(||). We consider a two-dimensional system with the 
gauge field A(r) described by the Hamiltonian H, Eqs.(0) 
and (|J). This Hamiltonian can be represented as the 
square of the Dirac operator n (we put below c = 1) 



W 

H = — = H - /i B cr 3 S 3 (r), 
zm 

n = a [p - eA(r)] . 



(3) 
(4) 



We assume that the distribution of the gauge field A(r) 
is Gaussian with the correlations 



2m 2 



S..-SW 

q 2 + K 2 p 2 F V y q 2 



(6) 



where 7 characterizes the strength of the disorder and pf 
is the Fermi momentum. The limit K = corresponds to 



R.A/ \ 

9 k (P, 



k±iS-ap + i: R ' A (p) : 



(10) 



where the self-energy S (p) should be found from the 
equation 

S^(p)=^m 2 X:y'^^(p-q) 



XCTi 



k±i6-aq + Y, R > A {q) 3 



a r (11) 



In Eqs.([10|) and (|l l|) , p and q are momenta. Using the 
usual assumption that disorder is weak one can obtain 
easily 



R.A/ \ 

9 k (P) 



1 



1 k + crp 

^ 2m x 2t 



(12) 
(13) 



2 



where r is the single-particle lifetime specified below. Us- 
ing Eq.(|^), we extract the Green function 

' (14) 



2m 2r 



This is the usual form of the Green functions. So, we 
conclude that the interaction of the magnetic spin does 
not change the form of the Green functions. 

At the same time, the expression of the single-particle 
lifetime r differs from that for spinless particles. For 
short-range disorder k 3> (Epr)^ 1 the solution of 
Eq.(|rH) leads to the following expression for the single 
particle lifetime t: 

7^Tg ■ .J, 3 (15) 



1 



•2tt 



sin 2 i + € 



2jE F 



(16) 



In contrast to, the corresponding result for the model 
without spintffl a factor cos 2 (6/2) in the integrand is ab- 
sent. However, in the limit of small k, integration over 
9 in Eq.(15) leads to the same result as for the spinless 
problem. 

In the limit of long-range disorder k <C (Epr)^ 1 , we 
obtain solving Eq.(|ll|) a more complicated expression 

1 2 1 E f t f°° Ae f 27r d9 . 2 9 



2tt 
1 



16EI 



AEt, 



lm 2 - 



2 

1/2 



16EI 



sin 



(17) 



(18) 



Again, the integrand in Eq.(|17]) differs from the corre- 
sponding integrand of the spinless problem by the ab- 
sence of cos 2 (9/2) in the integrand. This changes the fi- 
nal result and we obtain a single-particle lifetime r which 
has an additional factor 1 / v2 as compared to the corre- 
sponding result for the spinless problem. 

Due to the long range correlations of the random field 
the transport time Tt r entering the classical diffusion co- 
efficient D differs from the single particle time r. By 
considering the nonsingular corrections as explained in 
Ref. U one comes to the renormalization 

111 

(19) 



Ttr 



-(1) 



where can be obtained by inserting the additional 
factor cos 9 in the integrands in Eqs.([l5|) andr(|l7|). For 
the model without spin one has l/r tr = jEpQ and the 
diffusion constant is given by D = v F T tr /2. In the 
present case, one comes to the following expression for 
l/r tr [again, the factor cos 2 (9/2) is absent in the inte- 
grand] 



Ttr 



2jE F . 



d9 1 



-(1 



(20) 



The diffusion constant D is related to transport time Ttr 
in the usual way 



D 



(21) 



This result shows that by taking into account spin, we get 
a transport time which is one half of the usual one. This 
difference can be considered as a spin effect. Remarkably, 
the orbital and spin scattering give equal contributions 
to the resistivity determined by l/r tr . This effect is not 
trivial and is obtained only after integration over 9 in 
Eq.(p0[) and in the corresponding equation for the spin- 
less problem. It is specific for two dimensions. 



III. SUPERS YMMETRIC NONLINEAR a MODEL 

In order to consider interference effects and localization 
one should consider either more complicated diagrams or 
derive a nonlinear a model. Bpthlhc methods have been 
used for the spinless problem.E^ttl For the present prob- 
lem we want to use the supersymmetry technique and to 
derive a supermatrix a model. We use the notation and 
conventions of Rcf. ^ in the following calculation. 



As we have seen, we need the Green functions 



R,A 

9 k 



and g^' k A in order to calculate the Green functions G^' A . 
We define the following generating function to calculate 
the product of the Green's function as 



Z(k,Q) = J D^Dil)c- c , (22) 
£ = -i J d 2 r$(r) [H + eaSA(r)] rp(r), (23) 

where 



H = -iaaV — k+ —A. 



(24) 



k denotes the Fermi momentum and Co is related to the 
energy difference asw = WfU. The matrices S and a are 
defined as 



1 o 

-1 



CTi 
0~i 



(25) 
(26) 



Note that a are the (Pauli) matrices in the spin space, 
T3 in the time-reversal space and a is a matrix that 
reflects +IT and —II. The supervector rp contains 
32 components corresponding to fcrmion/boson, ad- 
vanced/retarded, time-reversal multiplication, spin de- 
grees of freedom, and +/— structure. 

The calculation is done in a similar way as that for the 
model without spin. After the ensemble averaging and 
the Hubbard-Stratonovitch transformation, the generat- 
ing function can be written as 



3 



(Z(k,u>)) = / DQexp 



J DQ cxp i str In ^iH Q 5(r - r') - 2m 2 ^ (r - r')aE,Q(r, r')aE 3 - j 
"IT^ / dWr%(r-r')str Q(ry)a£ i Q(r / ,r)aE J - 



(27) 



where Q is a 32 x 32 supermatrix and "str" denotes the 
supertrace. 

The saddle-point equation is solved by 



Q<°>(p) 



(28) 



which is essentially the same as Eq.@. The single- 
particle lifetime r has already been obtained above. 

Now let us turn to the discussion of the saddle-point 
manifold. For w = 0, the generating function is invariant 
under the transformation tp — > Tip, where the matrix T 
obeys the following constraints: 



TT = 1, 

= a: Ei 2- 



(29) 
(30) 



The first equation is the unitarity condition, whereas the 
second one specifies the structure in the spin space. Us- 
ing these conditions we obtain the saddle-point manifold. 
It can be parametrized as 



Q = TQMf, 



T = 



P = 



1-iP 
1 + iP 

q b 

B 



1/2 



B = 



(31) 
(32) 

(33) 



This block form of the matrix B reflects the structure in 
+/— space. This matrix satisfies [B, oE^a] = in ad- 
dition to the usual condition B = KBK. It means that 
and are proportional to the unit matrix in 

spin space and B^ ' and B^ ^ are proportional to E3. 
This result shows that if we consider the generating func- 
tion for the Dirac Hamiltonian n, only the mode 
is left gapless and we come to the usual unitary a model. 
The structure in +/— space makes the result more com- 
plicated. We have four gapless modes and they equally 
contribute to the physical quantities as we see below. 

Having specified the saddle-point manifold, we derive 
the free energy describing gapless modes. Again, this cal- 
culation resembles that for the spinless case. In addition 
to the usual massless mode -P(R), we need to take into 
account soft massive modes P(R, n) which depend on 
the direction vector of the Fermi momentum (n — p/|p|). 
The saddle-point manifold is parametrized as Eq. ( pl| ) and 
the matrix T is decomposed as 



T(R, n) = £7(R)y(R, n). 



(34) 



The supermatrices U and V parametrize the massless 
modes and the soft massive modes respectively. They 
are expressed as Eq.(|3^) and expanded in the superma- 
trix P. For the soft massive modes, the dependence of the 
supermatrices on n is treated in the Fourier-transformed 
space as 



F(R,n) =^P(R,m)e J 



(35) 



where 9 is the polar angle of the vector n. The mth life- 
time r'™) which is associated with each modes is defined 



as 



8£ - p , r (m) 



(1 + aSn)A 



= 2m 1 



i,3 



d 2 q 



(27T) 



2^ J (p-q)aE i QW(p)aE j e 



(o), 



(36) 



T (o) corresponds to the single-particle lifetime. We will 
see that these nonzero-harmonics modes renormalize the 
transport time and, as a result, one comes to Eq. (fl9|) . 
Details are presented in Ref. ||. Performing the gradient 
expansion and integrating out the soft massive modes we 
obtain the nonlinear a model 

F = l6~ / rf2rStr [^[ VQ(r)]2 + 2lujA Q( r )] . (37) 

where Q(r) — U(r)AU(r) and v = m/27r is the density 
of states. The diffusion constant D is given by Eq.([2l|) 
and the transport time Tt r is introduced in Eq. (po|) . 

The form of the free energy F is the same as for the 
model without spin but the dependence of the diffusion 
coefficient D on the correlations of the magnetic field is 
somewhat different. We note that an additional term de- 
rived in Ref. ^|can also be present in our model. However, 
this term appears at longer correlations of the magnetic 
field and we will not consider it in this paper. 

Due to the spin degrees of freedom, the 32 x 32 super- 
matrix Q has a more complicated structure than usual. 
However, we will show below, the renormalization prop- 
erties show this model falls into the usual unitary class. 
We show the results for the conductivity and the current 
correlation function. They can be calculated by using 
the following contraction rules for the perturbative cal- 
culations as 



4 



1 3 

(strXP(r)FP(r')) F = — H(r, r') ^[stra^Xstr a^Y - str Ac^Xstr Act,, Y 



16 

+str acr^Xstr aa^Y — str aAc^Xstr aAa^Y 
+strcr /i o'3Xstrcr AI (T3F — str Aa^a^Xstr Aa^a^Y 
—str acr M <T3Xstr aa^a^Y + str oA<7 M <73 Astr aAa^a^Y 

+ (X - T3 X,Y -> r 3 Y)], (38) 
1 3 

(strXP(r)strFP(r')) F = -n(r,r') ^ str (a^Xa^Y - Aa^XAa^Y + aa^Xaa^Y - aAa^XaAa^Y 

+<r fi <X3 Xa^aaY - Aa ^X Aa ^a^Y - aa^Xaa^Y + aAa^XaAa^Y) 
+(y-»-y), (39) 



where A and Y are arbitrary 32x32 matrices which 
commute with T3 and II (x, y) is the diffusion propagator 



Il(r, r ) = / - — 

V ' ' 2-kv J (2tt 



1 



a iq(r-r'; 



£>q 2 



(40) 



At first glance, the above expressions of the contraction 
rules look quite different from the usual ones. However, if 
we assume that A and Y do not have any spin structure, 
we obtain the standard contraction rules. Therefore, we 
will get the usual unitary results for spin-independent 
quantities. We show for instance the results of the con- 
ductivity and the current correlation function. The con- 
ductivity is defined as 



(41) 



where Py(k, to) is the Fourier-transform of the correla- 
tion function 

Py(r,r» - (tr^G^ , 2 (r,r')^,G« (r',r)\, 



(42) 

where 7r r = — iV r /m is the velocity operator. In the 
two-loop approximation we obtain 



(?ij(uj) = (7 6i. 



1 



2ir 2 v 2 



1 - 



d d q_ 



Z3q 2 — iuj 



(43) 



where <t = 2e 2 vD is the classical conductivity. This is 
the same expression as the one presented in Ref. || and is 
just the usual unitary result. 

The current correlation function which is defined as 

hi (r, r» = (Ji(r, E - w/2)Jj (v',E + w/2)) , (44) 



Ji(r, E) = — lim (tt* - 

47T «■'—"• 



X 



tr [G|(r,rO-G|(r,r')], 



(45) 



can be written in a Fourier-transformed form in the lead- 
ing order as 



ly(q,u;) = 



2e 2 



5ij - 



In (gLu,) , 



(46) 



where L u = (D/uj) 1 ^ 2 . This result is four times larger 
than the corresponding result of Ref. [l6|. This difference 
is due to the spin degrees of freedom. We see that all the 
expressions for the spinless quantities are just the results 
for the unitary ensemble. 

For the present model, we can introduce in addition 
spin correlation functions that could not exist for mod- 
els without spin. As an example, we define a correlation 
function which is a direct extention of the density-densitw 
correlation function widely used in theory of localization^ 



f dE 


f d2p \ 




(2tt) 2 [ 



(n(E) -n(E-w)) 

x (tra i G|(p)(7 J G^_ a ,(p-q)) 
+n{E) (traiGl+Jp + q)<r j G|(p) 

-tr^G^p^G^Jp-q)) 



(47) 



where n(E) is the Fermi distribution function and i,j = 
0,1,2,3 ((To = !)■ The component T o(q, oj) gives the 
conventional density-density correlation function. We re- 
fer to Appendix for details of the calculation and obtain 



iii(q,w) = 



2;/ 



_2al 



Dq 2 

2v 



for i — j — 0, 3, 
for i = j = 1, 2. 



(48) 



We find the usual result for the density correlation func- 
tion. Since we do not take any electron interaction into 
account, this can be considered a natural result. The 
same diffusive form appears for T33. This diffuson con- 
tribution comes from the mixing of the +/— space. 
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IV. CONCLUSION 
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In conclusion, we have discussed spin effects for the 
two-dimensional random magnetic held model. We sug- 
gested to take into account spin degrees of freedom writ- 
ing the Dirac Hamiltonian instead of the initial Hamil- 
tonian for electrons with spin, which is possible for the 
model of free electrons. The interaction of the magnetic 
field with spin is not smaller than its effect on the orbital 
motion. Moreover, in a homogeneous field the Zeeman 
splitting is equal to the distance between the Landau lev- 
els. 

We derived a nonlinear supermatrix a model that 
turned out to be a conventional a model with the uni- 
tary symmetry. However, the single particle lifetime and 
the diffusion coefficient differ from their values for the 
spinless particles. Remarkably, the transport time and, 
correspondingly, the diffusion coefficient are two times 
smaller than those for the spinless problem. This means 
that the orbital and spin scattering equally contribute to 
resistivity. The conductivity and the current-current cor- 
relation functions take the conventional form for the uni- 
tary class. The form of the spin-spin correlation function 
is the similar to that of the density-density correlation 
function. Our results obtained by the supersymmctry 
method can be well reproduced using the diagrammatic 
methods. 

Finally, we mention two possible extentions of our 
model. First, our analysis is restricted to the model with 
fixed g factor g = 2. The advantage to introduce the 
Dirac Hamiltonian was that the spin degrees of freedom 
is taken into account naturally. The square of the Dirac 
Hamiltonian without any parameter gives Eq.(^) with 
fixed g factor. On the other hand, it is well known that 
the g factor can take different value for realistic materi- 
als. Unfortunately, it is difficult to extend our approach 
directly to the case of a different g factor. However, we 
hope that our results shed some light on a possible be- 
havior of more realistic systems. 

Second interesting problem is the case of a nonzero av- 
erage magnetic field. In this case, the Zeeman splitting 
causes the reformation of the Landau levels and nontriv- 
ial results may be expected. For spinful electrons, in 
addition to the perpendicular magnetic field, the parallel 
field is also relevant even for the two-dimensional sys- 
tem. In this paper, we have discussed a two-dimensional 
Hamiltonian. The Dirac Hamiltonian does not include 
the Pauli matrix 03 and only the perpendicular mag- 
netic field enters the square of the Dirac Hamiltonian 
[Eq.@]. If we consider the three-dimensional Hamilto- 
nian, we can include the magnetic field for all directions. 
Even if the diffusion process is two-dimensional, we may 
need to treat the one-electron states as three-dimensional 
in that case. It will be a subject of a future work. 
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APPENDIX A: CALCULATION OF THE SPIN 
CORRELATION FUNCTION 



We consider the function 

(tr^Gg^rV.G^r',!-)), 



(Al) 



to calculate the spin correlation function. For the prod- 
uct of the Green functions, we express it in the functional 
integral form as 



G^>,r')G^(r',r) 



= / D(^)[^+(r)^(r)+V^(r)^(r)] 

V F J 

x[^+(r')^( r ') +^(r')V^(r')]e" £ , (A2) 

where a, (3, 7, and S are spin indices, 1 and 2 denote the 
advanced and retarded channel, + and — the +/— space, 
and n and v are other indices. No summation is implied 
for (i and v. After introducing the auxiliary field Q and 
taking the contraction of ip in two possible ways, we get 



1 



(cr l ) Q/ 3(er,-)7<5 str C+gp a (r, v)C-g Sl {r' , r' 



~~T\ (e^WOj^str G+g^ 7 (r,r') str C-g Sa (r',r) 



where 



C - k 1 ~ T3 1 ± A 
± 2 2 2 



(A3) 



(A4) 



At the leading order, the Green's function g is given in 
the Fourier-transformed space in terms of the relative co- 
ordinates as 



s(R,p) 



aSp-e + ^(e + a£p)Q(R)' 



(A5) 



and Q(R) is expanded in terms of P(R). 

Substituting g to the first term of Eq.(A3), we have 
leading contribution 



(cTj) 7 «5Str \C + P l3a {Y)C-Ps 1 {Y l )\ ) 



The matrix P has the structure 
p(++) p(+- 



P = 



p(-+) pi- 



rn 



(A7) 



G 



in +/— space. p(++) and p( __ ) are proportional to ijq 
in spin space and P^ ' and p( ^ are proportional to 
03. Using this fact and the contraction rule, we get 

2tt V 6 i0 6 j0 (sbc C+P {++) (r)CUP (++) (r') 

+str C+P { — ] (t)C-P { — ) (r')\ 
+5 i3 ^3(strC r + p(+-)(r)C_p(-+)(r / ) 

+strC+P ( - +) (r)C-P (+ - ) (r')) 
d 2 q 



47rr/ (5jo^o + ^3^3) 



(2tt) 2 L>q 2 



iq(r-r') 



(A8) 



The contribution comes from the second term of Eq.(A2) 
is small and is neglected. 

We use this result for the first term of Eq.([47j). Sub- 
stituting the Green's function ( |l4| ) to the second term 
of Eq. (|47|) , we find 2v5ij. Combining these results, we 
finally obtain Eq.(p|). 
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